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Abstract: We prove two limit relations between Sclirodinger operators 
perturbed by measures. First, weak convergence of finite real- valued 
Radon measures /i„ — ^ m implies that the operators —A + e^A^ + /i„ 
in LF'{W^,dx) converge to —A + e^A^ + m in the norm resolvent sense, 
provided d < ?>. Second, for a large family, including the Kato class, of 
real- valued Radon measures m, the operators — A + e^A^ + m tend to the 
operator —A -|- m in the norm resolvent sense as e tends to zero. Explicit 
upper bounds for the rates of convergences are derived. Since one can 
choose point measures /i„ with mass at only finitely many points, a com- 
bination of both convergence results leads to an efficient method for the 
numerical computation of the eigenvalues in the discrete spectrum and 
corresponding eigenfunctions of Schrodinger operators. The approxima- 
tion is illustrated by numerical calculations of eigenvalues for one simple 
example of measure m. 



I Introduction 

In this paper we are going to analyze convergence of Schrodinger operators perturbed 
by measures. It is known that weak convergence of potentials implies norm-resolvent 
convergence of the corresponding one- dimensional Schrodinger operators. This re- 
sult from 1^ may be interesting for several reasons. For instance every finite real- 
valued Radon measure on M is the weak limit of a sequence of point measures with 
mass at only finitely many points. There exist efficient numerical methods for the 
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computation of the eigenvalues and corresponding eigenf unctions of one-dimensional 
Schrodinger operators with a potential supported by a finite set; actually the effort 
for the computation grows at most linearly with the number of points of the support 
P]. Since the norm resolvent convergence implies convergence of the eigenvalues in 
the discrete spectra and corresponding eigenspaces, we get an efficient method for 
the numerical calculation of the discrete spectra of one-dimensional Schrodinger op- 
erators. Norm resolvent convergence has also other important consequences: locally 
uniform convergence of the associated unitary groups and semigroups, convergence of 
the spectral projectors (which implies the mentioned results on the discrete spectra) 
etc. 

Let us also mention a completely different motivation for studying convergence of 
operators with point potentials. In quantum mechanics neutron scattering is often 
described via so called zero-range Hamiltonians (the monograph pP is an excellent 
standard reference to this research area). In a wide variety of models the positions of 
the neutrons are described via a family (Xj)"^^ of independent random variables with 
joint distribution fi. Usually the number n of neutrons is large and one is interested 
in the limit when n tends to infinity and the strengths of the single size potentials 
tend to zero. In the one- dimensional case this motivates to investigate the limits of 
operators of the form 

a 7^ being a real constant and {Q, JF, P) a probability space. By the theorem of 
Glivenko-Cantelli, for P-almost all u E Q the sequence (^ Yl]=i^Xj{oj))nm converges 
to the measure a/j, weakly. By the mentioned result from 6j, this implies that 



in the norm resolvent sense P-a.s. 

It is the purpose of the present note to derive analogous results in the two- 
and three-dimensional case. It was shown in [6] and [8^ that one can approximate 
Schrodinger operators perturbed by suitable measures by point potential Hamilto- 
nian. However, the convergence there was in the strong resolvent sense, which is of 
course a weaker result than the norm resolvent convergence. 

If the dimension is higher than one, then it seems to be impossible to work directly 
with operators of the form —A + /i, n being a point measure. In fact, while the 
operators — ^ + Yl^=i(^j^xj can be defined in dimension one via Kato's quadratic 
form method as the unique lower semibounded self-adjoint operator associated to the 
energy form 

D{S) := H\R), 

/n 
|/'(x)prfx + 5^a,|/(x,)p, feD{S), 
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/ being the unique continuous representative of / G if^(R), in higher dimension 
(i > 1, the quadratic form 

D(£) := {/ G H^{W^) : f has a continuous representative /}, 



S{f,f) ■■ = 




is not lower semibounded and closable if at least one coefficient aj is different from 
zero. 

The strategy to overcome the mentioned problem in higher dimensions is based 
on two simple observations: 

1. The lower semibounded self-adjoint operator + /i can be defined via Kato's 
quadratic form method for every real-valued finite Radon measure /i on M'^ (if d G 
{1,2,3}), including point measures. 

2. —A + 5^A^ — > —A in the norm resolvent sense, as e > tends to zero. 

We show the convergence claim in two steps. In section |n] we shall prove that the 
sequence (—A + e^A"^ + /i„)„gN converges to —A + e^A"^ + m in the norm resolvent 
sense provided d < 3, e > and the finite real- valued Radon measures fin on 
converge to the finite real-valued Radon measure m weakly. Then, for a large class 
of measures m we shall prove that 

-A + £:^A^ + m — > -A + m 

in the norm resolvent sense as e tends to zero, cf. section IIIIl Actually, we will not 
only prove convergence but also give explicit error estimates. 

As approximating measures /i„ we can, in particular, choose point measures with 
mass at only finitely many points. In section HVl we will present formulae which make 
it possible to calculate the eigenvalues and corresponding eigenspaces of operators 
perturbed by a finite number point measures. Then similarly to ^ chapter II. 2], 
the spectral problem means to solve an implicit equation and the effort for these 
computations grows at most as 0{n'^). 

Putting both convergence results from sections |n] and IIIIl and formulae from sec- 
tion llVl together. we get an efficient method to calculate the eigenvalues in the discrete 
spectrum and corresponding eigenspaces of Schrodinger operators —A + m numer- 
ically. We apply the approximation to the simple two-dimensional example, where 
measure m is negative and supported by a circle. 

Our method does not only cover the case when m is absolutely continuous w.r.t. 
the {d — l)-dimensional volume measure of a manifold with codimension one but a 
fairly large class of measures m containing the set of all finite real-valued measures 
belonging to the Kato class. In particular, the absolutely continuous case dm = Vdx 
where —A + m = —A + V is a. regular Schrodinger operator is contained in our 
approach. We refer to [TU] for related convergence results in the regular case. 

Notation and auxiliary results: Let yU be a real-valued Radon measure on M'^. 
By the Hahn- Jordan theorem, there exist unique positive Radon measures /j,^ on 
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such that 



^ = ^+ and /i+(R'^ \B) = = ^.-{B) 
for some suitably chosen Borel set B. We put 

II II := /i+(R'^) + and |/i| := /i+ + /i". 

If /i is finite, then we define its Fourier transform fi as 

A(p):^(2.)-/^/e-,(..). peR^. 

Similarly, / also denotes the Fourier transform of / G L'^{dx) := L^(M'^, (ix), rfx being 
the Lebesgue measure. 

For s > we denote the Sobolev space of order s by if^(]R'^), i.e. 

H'iR') := S^feL^dx): J (1 + py\fip)\'dp < 

WfllHs := (^Jil+p'r\fip)\'dpy , feH^iR"). 



We shall use occasionally the abbreviations := L^(]R'^,/i) and H'^ := if^(M'^). 

II T \\ni,n2 denotes the operator norm of T as an operator from Tii to 7^2 and 
II T ||'H-=II T \\n,H- II / \\h and (/, represent the norm and the scalar product in 
the Hilbert Ti, respectively. If the reference to a measure is missing, then we tacitly 
refer to the Lebesgue measure dx. For instance "integrable" means "integrable w.r.t. 
dx" if not stated otherwise, || T ||, {f,h) and || / || denote the operator norm of T, 
scalar product and norm in the Hilbert space L'^{dx), respectively. We denote by 
C^(R°') the space of smooth functions with compact support. 

For arbitrary e > {e = will be admitted only in section IIII|) let Se be the 
nonnegative closed quadratic form in the Hilbert space L'^{dx) associated to the 
nonnegative self-adjoint operator —A + e^A^ in L'^{dx). Obviously we have 

SeUJ) = £'(A/,A/) + (V/,V/)>£2(A/,A/), feD{Ss), 

for every e > 0. Note that for 6 = the form domain is H^{R'^) and So is the classical 
Dirichlet form. For any a > we put 

£Uf^ h) := h) + h), f,he D{£,). 
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II Operator norm convergence 

Throughout this section let d < 3 and /i be a finite real-valued Radon measure on 
M'^. Then, by Sobolev's embedding theorem, for every s > 3/2, and, in particular, 
for s = 2, every / G if*(M'^) has a unique continuous representative / and 

II / |U:= sup{|/(x)| : X G M'^} < c, II / \\hs, f G H'(R''), (1) 

for some finite constant Cg. Note that < 1 if s = 2. It follows that for every e > 
and every r] > there exists an a = a{e, rj) < oo such that 

\\f\\lo<v£eifJ) + aifJ), feH'{R% (2) 

Since fi is finite, for arbitrary 6,7] > and some finite a we get 

I J |/|>| <v\\f^\\ £eU, /) + « II II (/, /), / e H\R'). (3) 

We put 

By and the KLMN-theorem, is a lower semibounded closed quadratic form in 
L'^{dx). We denote the lower semibounded self-adjoint operator in LF'{dx) associated 
to £^ by -A + e^A^ + /i. 

Our main tool to prove convergence results will be a Krein-like formula which 
expresses the resolvent (—A + e^A^ + /i + a)~^ by means of the resolvent 



G,., := (-A + £2A2 + a^-l 



The operator Ge,a has the integral kernel ge,a{x — y) with Fourier transform 

1 

^Zp4 _|_ _j_ Q, ' 



For every e > and a > 0, the function g£.a{x) is continuous on \ {0} and if d = 1 
or if d < 3 and £ > it is continuous on whole W^. Moreover, it is radially symmetric. 
Finally, (7o,a is the Green function of the free Laplacian in R"' and it is nonnegative. 
By the dominated convergence theorem, 

II 9e,a \\]l2= / I o 4 I 2 I \2 ' 0' ^ l"l ' ^ 

which, by Sobolev's inequality, implies that 

II 9e,o, lloo — ' 0, as |a| — > oo. (5) 
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The fact that Qe^a is the Green function of —A + e^A^ means that 

j ge,a{x — y){—A + e'^A'^ + a)h{y)dy = h{x) dx-a.e. 

for all h G D{—A + e^A'^) = if^(M'^). The equation above does not only hold almost 
everywhere w.r.t. the Lebesgue measure dx but even pointwise everywhere, as the 
following lemma states. 

LEMMA 1 Let Green function ge,a o-nd operator — A -\- e'^ A'^ + a be defined 
as above. Then one has 

j ge,a{x - y)i-A + £^A^ + a)h{y)dy = h{x), x eW^ (6) 
for all h G if^(M^). 

Proof: In fact, we have only to show that the integral on the left hand side is 
a continuous function of x G M'^. We choose any sequence (/n)neN of continuous 
functions with compact support converging to (—A + e^A^ + a)h in L'^{dx). By (jlj), 
Qe^a € H'^{M.'^) C L'^{dx), therefore we can write 

geA^-y)i-^ + ^^^^ + (^)h{y)dy= lim geA^-y)fn{y)dy, xgR'^. 
J — "^J 

Obviously the mapping x J gs^x — y) fn{y)dy , M"' — > C, is the unique continuous 

representative Ge^afn of Gir^afn for every n G N. Since G^^a is a bounded operator 
from L'^{dx) to H'^{R'^) (even to if^(M'^)), the sequence {Ge,afn)nm converges in 
to Ge,a{-A + s^A'^ + a)h = h. By Sobolev's inequality JH), this implies 

that the sequence (Ge,a/n)neN of the unique continuous representatives converges to 
a continuous function uniformly. By the last equality, x J g^^x — y){—A + e'^A'^ + 
a)h{y)dy, M*^ — > C, is this continuous uniform limit and we have proved ©• n 
We introduce following integral operator 

G^Jix) := j gU^ - y)f{y)f^{dy) rfx-a.e., / G H\R''). 

We can prove several estimates of its operator norm. 

LEMMA 2 The operator G'^ ^ is bounded on if^(M'^) and its operator norm 
II G^Q, 11^:^2 decays with a — ^ oo. The operator is bounded also w.r.t. other operator 
norms, in particular there are finite real numbers Cj, i = 1, 2, 3 such that 

Wg^jWh-^ < ci(«)||/iu 

IIG.VIU^ < C2(a) II / lU.(H) feH^R") 
II G^,af ||L2(|/i|) < C3(a) II / ||l2(|m|) 
and all three numbers Ci vanish in the limit a — > oo. 
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Proof: Using Sobolev's inequality we have for arbitrary / G if^(M"') 

\Mp)\' < (2vr)-'^ II / IILII /i r< {27r)-' \\ f \\U\ f , P ^ 
Then the convolution theorem yields 

WG^wj,. = J\{i+pr\u,a*h){p)\'dp 



- I 1-2^4 ^^2 ^^,12 II / ll^ll ^ II' 



|^2p4 _|_ p2 _j_ ^ 



< / I V, i. ^P ||/|||.||Mf<oo, 



Therefore G^^, is an everywhere defined bounded operator on i^^(R'^) and we get an 
upper bound for the norm 

and the expression on the r.h.s. is also the uniform upper bound Ci. 
To determine the remaining upper bounds C2 and C3, we can write 

G':j{x)\'dx 

QsA^ - y)fiy)fJ'^idy) - / g^^x - y)f{y)fi''{dy)\'^dx 



< J ^ J 9sA^-y)fiy)l^^idy)\'^dx + 2 J \ J g^A^ - y)f{y)n {dy)\^dx 

< 2 y y |(7,,,(x-y)|V(c?Z/) j \f{y)\^ii-^{dy)dx 
+ '^J J \9eA^-y)\'f^'idy) J \fiy)\'fi-idy)dx 

< 2 J \gU^)\'dx II II J \f{y)\'\f,\{dy), f G H\R'). (8) 
In a similar way we arrive at 

\GCf{x)\'\f,\{dx) < 2 II g,,^ Will ^ f /" \f{y)\'\fi\{dy). 



Finally, from (jH) and (0) one concludes that all the upper bounds of the operator 
norms tend to zero in the limit a — > 00. □ 
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General results of jS] (cf. also section UTTl below) provide, in particular, an explicit 
formula for the resolvent of the operator —A + e^A^ + fi. In this resolvent formula 
there occur operators acting in different Hilbert spaces. This is inconvenient when we 
investigate the convergence of sequences of such operators and we shall use a slightly 
different resolvent formula: 

(-A + e^A^ + + a)-' = G,,, - G^JI + G^J^'G,,^. (9) 

For the sake of completeness we present the proof of the above Krein's formula in the 
appendix. According to lemma |2l we can choose a > such that || G^,^ \\h^,h'^< 1- 
Then the operator / + is invertible and its inverse is everywhere defined on 
if^(R'^) and bounded; here / denotes the identity on //^(M*^). By 0, we can choose 
a > such that, in addition, 

S^JfJ):=S^{fJ) + a{fJ)>{fJ), feD{S^). (10) 

We are now prepared for the proof of the main theorem of this section: 

THEOREM 3 Let m and fin, n G N, be finite real-valued Radon measures 
on M.'^. Suppose that the sequence (/in)ngN converges to m weakly and sup„gj^ || Hn \\< 
oo. Let e,a > and d G {1, 2, 3}. Then the operators —A + e^A^ + yU„ converge to 
—A + e^/S? + m in the norm resolvent sense. 

Proof: Let e > be arbitrary. We choose < c < 1 and a > such that 

II f^nf [ I .l^^fi' I, dp <c\ ne M, (11) 
J \e^p^ + p^ + 

and 

ll"^f / i2^!l^Pl ,. dp<c\ (12) 
J \e^p^ + p^ + a]-^ 

According to Q, we can choose a > such that, in addition, 

Sta{fJ)>{fJ), feH'iR"), neK (13) 

Since (/i„)„eN converges to m weakly, (fT^ also holds when we replace fin by m. By 
Lemma 121 in particular estimate (|7j), inequalities pi|) and (fT^ yield 

\\G^^-^\\h2,h^ < c, nGN, 

\\GT^J\h2,h^ < c, (14) 
||g::,/||h^ < c||/|U, feH\R'). 

Hence the resolvent formula is valid both for /i = m and for /i = n G N. By 
Lemma 121 we can choose a sufficiently large so that also 

|G™Xx)pc/x < |/i|2(i|m| and j \G^{x)\^\m\{dx) < j \h\^d\m\ (15) 



for every h e H'^(R'^). 

For notational brevity we put 

^7o:=^7o,i, g-=ge,a, G:=G,,^, G^" := G^,; and := 
With this notation we have 

(-A + e^A^ + ^„ + - (-A + e^A^ + m + 

= (G" - G^")[/ + G'^J-iG + (G'^" - G™)[/ + G'"]-i(G'^" - G™)[/ + G^-J'^G 
+G"[/ + G'"]-i(G^" - G'")[/ + G'^"]-iG. 

Since G is a bounded operator from L^{dx) to if^(]R^) we have only to show that 

II G'" - G^" Wn^L^d.)-^ as n ^ oo, (16) 

II G^[I + GH^'IG™ - G^") ||^,2,^2(,,)— . 0, as n ^ oo. (17) 
We introduce 

z/„ := m - /i„ 
i^nx(c^l/) := gix - y) Vnidy), X e M^, n G N. 
As (i < 3, the function 

is continuous and bounded for every / G L?'{dx)\ this well known fact can be proved 
in the same way as (P). Since the function g is bounded and go is nonnegative it 
follows that 



9{x-y)\ / l5'o(z/ - a)| |(-A + l)/i(a)| rfai/^(c/?/) 



< oo 



for all a; G M"' and h G //^(M"'). Hence by Fubini's theorem, the function k,^^^ : R'^ 
R, defined by 

Jgoiy-a)g{x-y)unidy), if defined, 
0, otherwise. 



is Borel measurable, the integral on the right hand side is defined and finite for almost 
all a G M*^ (almost all w.r.t. the Lebesgue measure) and 



\{G^-h){x)\' = I / g{x - y) h{y) ur,{dy)\' 

9{x-y) / go{y - a){-A + l)h{a)daUn{dy)\'^ 



< J \K„Aa)\^da- J \{-A + l)h{a)\^da 

J \K^^{a)\^da, he H''{R''),neN. (18) 



< 2\\ h \\l 



Thus in order to prove ()16|) we have only to show that 

\k„„^{a)\'^dadx — >0 as n — ^ oo. (19) 




We have 

/ J \ku„Aa)\'^dadx = i2nY J J \go{p)\'^\uZ{p)\'^dpdx 
= 1! jYT^ J e'Pygix - y) u^) I e~'^^g{x - dpdx. (20) 

Since |1 and g are integrable w.r.t. the Lebesgue measure, g is bounded and 

the Radon measures z/„ are finite, we can change the order of integration. Let us 
rewrite as 



/(y, z) h{y, z) ® Vn{dydz). 



The function 



/(y, z) := j e'Pye-'^^j^-^^ dp, y, z e R", 

is bounded and continuous. It follows from the fact that it is (up to multiplication 
by (2'7r)'^/^) the inverse Fourier transform of the integrable function |1 at the 

point z — y. 

Also the function 

h{y^ ^) := / 9{x - y) g{x - z) dx 



is bounded and continuous for y,z & M.'^. This can be shown using following observa- 
tion. Let ?/ G M'^ and K be any compact neighborhood of y. Since \x\^gi;^a{x) — > 
for every j G as |a;| — > oo, there exists a constant a < oo such that 

\g{x-y) g{x- z)\ < a \\ g Woo dist{x,K)-'^, xeR'^\K, z eR'^, y e K. 

By Stone- Weierstrass theorem, the set of functions of the form X^jLi fj{^)9'j{y)i 
G N, where fj,gj are bounded and continuous, is dense in the space of bounded 
continuous functions w.r.t. the supremum norm. Since the measures tend to zero 
weakly and sup„gp^ || z/„ ||< oo, this implies that the product measures Vn ® I'n tend 
to zero weakly, too. Hence by (pUj) . we have proved (fT^ and therefore also (jTBI). 
It only remains to prove pTj) . For this purpose we first note that 




\kv„A'^)\ da\m\{dx) — >0 as n — > oo. 
This can be shown by mimicking the proof of (fT^ . By (jl8|) . it follows that 
j |(G'^"/i)(x)|2 \m\{dx) < 2cn II h 11^2, h G H^{W^). 
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Thus, in order to prove (fT7j) . we have only to show that there exists a finite constant 
C such that 

II G^{I + G^-^h |U2(,,)< c(^j \h\H\m\^ ' , he H^R"). (21) 
Using the estimates (fT^ . we have 



oo 



G™(J + G"")-^ = - ^(-G"")^'. (22) 



According to (jlSp . 



y |(G'™)i/i|2rf|m|J <c-c' fy |/i|'c?|m| 

for every j G N and hence 

oo oo . X 1/2 r f r ~ \^/^ 

II 5^(-G"")^7i |U^(d.)< 5^c^' \h\''d\m\^ ^T^c[j '^''^'"'IJ • 

By fl22p . this imphes ()2H) and the proof of the theorem is complete. □ 
REMARK 4 We have shown that 

II (-A + e^A^ + /i„ + a)-i - (-A + e^A^ + m + a)'^ f 

+^2 / / I / 5'o,i(l/ - a)9eA^ -y){m- fin){dy)\'^da \m\{dx) 

for some finite constants Cj = Cj{e,a), j = 1,2, which can be computed with the 
aid of the proof of theorem El Thus the proof provides explicit upper bounds for the 
error one makes when one replaces the operator —A + A^ + m by —A + A^ + /i„. 

REMARK 5 The essential spectrum of —A + A^ + m remains the same 
for any finite real- valued Radon measure m on R"' 

aess{-A + e^A^ + m) = (Te,,(-A + e^A^) = |0, oo). 

By Sobolev's inequality and [U Lemma 19], the mapping f ^ f from H'^{M.^) to 
L^(|m|) is compact. Therefore using estimate (jSj), one may conclude that G^^, is 
compact if regarded as an operator from if^(]R'^) to LF'{dx). According to the resolvent 
formula Q, this implies that the resolvent difference G^^ — G^.a is compact and hence 
the corresponding essential spectra coincide. 
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Ill Dependence on the coupling constant 



In this section we are going to prove that 

- A + er^A^ + m — > - A + m as e j 0, (23) 

in the norm resolvent sense. Here m denotes a real- valued Radon measure on R'' and 
we assume, in addition, that for every r] > there exists a < cx) such that 



\f?d\m\<i^(^j \Wf\'dx + ^J Ifl'dx 



f G C^{W). (24) 



Note that we neither require that m is finite nor that d < 3. On the other hand, the 
condition (j^^ implies that m{B) = for every Borel set B with classical capacity 
zero and, for instance, it is excluded that m is a point measure if > 1. 

The inequality holds, in particular, provided m belongs to the Kato class, 

i.e. 

sup |m|([n, n + 1]) < oo, d=l, 
lim sup / I log(|a; — y|)| |m|((i?/) = 0, d = 2, 



f 1 

lim sup / \m\{dy) = 0, d = 3 

xm^ J B{x,e) F - 2/1 



with B{x,e) denoting the ball of radius e centered at x (cf. ^T], Theorem 3.1). We 
refer to [3 chapter 1.2], for additional examples of measures satisfying (pH) . 

In general, the elements / in the form domain of —A do not possess a continuous 
representative /. Therefore we shall give a definition of different from the one in 
sectional so that it works for all e > 0. Of course, both definitions are equivalent in 
the special case of positive e. 

Since the space C^(]R'^) of smooth functions with compact support is dense in the 
Sobolev space H^{M.'^), there exists a unique bounded linear mapping '■ -ff^(M°') — 
L^(|m|) satisfying 



Jmf — fi / G CcT 



(strictly speaking maps the cix-equivalence class of the continuous function / G 
C^(R°') to the I m|-equi valence class of /). We put 

D{8T) := D{£,), 
STUJ) ■■= UfJ) + {AM,Jmf)LH\m\). feD{S^ 
where 0(8^) = H\M.'^) for e = 0, D{Se) = H^{R'^) otherwise and 

A^h{x) := i ^^^]\ "^^fd. D heLWm\), 
^ ^ 1 -h{x), xeW^\B, ^' 
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with B being any Borel set such that m"'"(]R'^ \ B) = = m^{B). By ()24|) and the 
KLMN-theorem, the quadratic form E"^ in L?'[dx) is lower semibounded and closed 
and 

^,':^,(/,/)>0, feD{Sr). 

Again, — A+e^A^+m denotes the lower semibounded self-adjoint operator associated 
to and we put 

K^^^ := i-A + E^A^ + m + a)-' 

provided the inverse operator exists. G^^a is defined the same way as in section HTl 

One key for the proof of the convergence result (|23|) is the observation that one 
can decompose 

. , , c(e) c(e) 
9eAP) 



whenever c{e) is defined. The coefficients —a{e) and —P{e) are the roots of the 
polynomial e'^x'^ + x + a; a simple calculation yields 

c{e) := ^ — ^ 1, as e i 0, 

2a 

a{e) := = — > a, a.s e [ 0, (25) 



1 + VI - 4:6' 



a 



(3{e) := ^ > oo, as £ J, 0. 



1 + VI - 4£2a 
Using the parameters introduced above, we arrive at 



Gs,a = c{e)Go,a{e) " c(£:)G'o,/3(e) • (26) 

In the proof of the convergence result (^Hj) we will use again a Krein-like resolvent 
formula, this time using the one from cf. ()28|) below. First we need some prepara- 
tion. Let a > and e > 0. We introduce the operator Jm,e,a from the Hilbert space 
{D{Ss),Se^a) to L^(|m|) as follows: 

Jm,£,af ■ Jmfi f ^ ^(^Jm,e,a) ■ 

By the operator norm of Jm,e,a is less than or equal to t] provided a > Thus 
we can choose ao > and c < 1 such that 

II Jm,£,a \\ (D {£e) ,£e .a) (Iml) < Vc, a> ao- (27) 

Due to ()27p. the hypothesis of Theorem 3 in Pj is satisfied and the theorem implies 
that —a belongs to the resolvent set of —A + e^A^ + m and 

a > tto- (28) 
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In fact, we can write 

J*m,e,a' = (JruGe^a'T, «' > 0, (29) 

since we have 

for every h G L'^{dx), s > and a' > 0. 

THEOREM 6 Let m be a real-valued Radon measure on M.'^ satisfying |^[). 
Then the operators —A + e^A^ + m converge to — A + m in the norm resolvent sense 
as e I 0. 

Proof: Both resolvents are written by means of Krein's formula ()28|). so we can 
compare the first and second terms separately. To see that || Ge^a — Go^a Wi'^idx) 
vanishes in the limit £ | is simple. It is enough to use the first resolvent formula, 

Go,a{e) — Go^a = {a — a{e))Go^aGo^a(e) (30) 

and the fact that 

for some continuous function k vanishing at infinity (actually, k{x) = for x <2 
and k{x) = l/(4(x — 1)) for x > 2). Then the decomposition of Ge,a and the 
asymptotic behavior of a{e),P{e) and c(e) finish the argument. 

The proof that also the difference of second terms in Krein's formula tend to zero 
as £ ^ can be reduced into two tasks 

II JmGe,a — JmGo^a \\L'2{dx),L2{\m\) ^ aS £ | 0, 

II (1 + JmJm,e,a^m) ^ ~ (1 + JmJmfl,a^'rn) ^ lU^dml) ^ aS £ | 0. 

The argument for the first line is similar to the one we have presented above for 
Ge,a — Go^a, wc ouly havc to add that, by hypothesis (j21|), it follows that 

II JmGo, «' lli2(dx),L2(H)< max(l,/5i)/c(a'), a' > 0, (31) 

where function k{a') is defined as above. 

To show the second line we choose any a > ao, then from (P7j) we get 

II (1 + JmJm,e,a^rn)~'^ |U2(|mj)< e>0. 

By the second resolvent identity 

(1 + A)-' - (1 + B)-' = (1 + A)-\B -A)il + B)-\ 
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it is sufficient to prove that 



II ■JmJm,£,a ~ ■JmJm,0,a \\L^{\m\) ^ aS £ j 0. (32) 

From (j2Sl) and ^ follows that 

JrnJm,e,a ~ '^i^) J m{J mG Q ,a{e))* — c(£:) Jm(<^mG'o,/3(e) )* ; 

note that c[e) is real for sufficiently small e. Using this expression and ()30|) and ()29|) . 

we get 

II T T* T T* II 

II '^rn'Jm,£,a ■^ni'^m,0,a ll-L2(|m|) 
— II (c(^) ~ l)^m(<^mG'o,Q(e))* ||l2(|„|-) + || J m{J rriG Q ,a{e))* — Jm{JmGo^a)* ||L2(|m|) 

+ II c{e)Jm{JmGo,l3{e))* lU^dml) 
= II — ^)Jrn,0,a(e)Jmfi,a{e) IU^(|m|) + II ^ « ) -'m.G'ca ( JmG'o,a(£) )* ||L2(|m|) 

+ II cie)Jm,Q,f3ie)J^^o,l3{e) \\L^{\m\), € > 0. 

According to (j2D), the mapping || Jm,o,aJmOa \\L^{\m\) is locally bounded for a G 
(0, oo) and tends to zero as a tends to infinity. Since a{e) — > a, c{e) — > 1 and 
P{e) — > oo as £ I 0, this implies, in conjunction with (jSD), that (jH^ holds. □ 

REMARK 7 By the proof above, || - || is upper bounded by an 
expression of the form c ■ (e^ + ri{m, e)) where the finite constant c can be extracted 
from the proof and r]{m, e) has to be chosen (and can be chosen) such that ()24|) holds 
with r] and P replaced by ri{m,e) and P{e), respectively. 

IV Eigenvalues and eigenspaces of the approxi- 
mating operators 

Throughout this section let d < 3 and let m be a finite real-valued Radon measure 
satisfying ()24|) (e.g., let m be from the Kato class). By the two preceding convergence 
results, we can approximate the operator —A + m in L^(M'^, dx) by operators of the 
form —A + e^A'^ + /i, where e > and /x is a point measure with mass at only 
finitely many points. Since the convergence is in norm resolvent sense, we can thus 
approximate the negative eigenvalues and corresponding eigenspaces of the former 
operator by the corresponding eigenvalues and eigenfunctions of the latter one. Note 
that we know from remark El and |5l Theorem 3.1] that the essential spectra coincide. 

The following theorem shows how to compute the eigenvalues and corresponding 
eigenspaces of the approximating operators. 

THEOREM 8 Let d < 3 and e > 0. Let fi = ^^^^^ Cj6^., where N eN, 
Xi, . . . ,xn are N distinct points in M'^ and Ci, . . . , cat are real numbers different from 
zero. Then the following holds: 
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a) The real number —a < is an eigenvalue of —A + e^A^ + if and only if 

det ( — + geA^j - ^k)j = 0. 

y'^k J i<j,k<N 

b ) For every eigenvalue —a < the corresponding eigenfunctions have the following 
form 

^ hkgeA- - ^k), {hk)l<k<N e ker I ^ + gsA^j - ^k) I 

k=l ' ^<j,k<N 

Proof: Since D^E^J) = H'^{R'^), the mapping can be understood as 

J J := / |/i|-a.e., / G H^M.''). 

By (0), / deA' ~ y)fiy)d-y is the unique continuous representative of G^^af ■ Hence 
J^Gs^a is the integral operator from L'^{dx) to L^(|/i|) with kernel g^A^ ~ v) 
inverse operator {J^G^A* is the integral operator from L^(|/i|) to L'^{dx) with the 
same kernel. Thus we get 

N 

J,,{J,j,GeA*^i^K^j) = ^CkgsA^j - Xk)h{xk), l<j<N, (33) 

k=l 

for every h G L^(|yu|). 

Due to Krein's formula (PH|) . —a < belongs to the resolvent set of (— A+e^A^+yu) 
provided 1 + J^{J^GsA*^ti i^ bijective. Since L^(|/i|) is finite dimensional and we 
have expression (j^Hj) . that is true if and only if 

A(a) := det{6jk + CkgeA^j - Xk))i<j,k<N 7^ 0, 

with Sj^k being the Kronecker delta. As geA^) is a real analytic function of a G (0, oo) 
for every x G M'', the function A(a) is also real analytic on (0, oo). By it is different 
from zero for all sufficiently large a. Thus the set of zeros on (0, oo) of this function 
is discrete. 

Since Jfj,Ge,a is surjective and {J^GsA*^fj. injective, the resolvent formula (PHj) 
implies that any ao > satisfying A(ao) = is a pole of R^a- Thus we have proved 
that — ao is an eigenvalue of —A + e^A^ + yU if and only if A(a;o) = 0. Finally, the 
expression 

det{6jk + CkgeA^j - Xk))i<j,k<N = nf^^Cfc ■ det{Sjk/ck + geA^j - Xk))i<j,k<N 

implies the assertion a). 

By the preceding considerations and [SI Lemma 1], 

h ^ {J^G,A*\h 

is a linear bijective mapping from ker(l + J^G'g^Q,)*^^) onto ker(— A+e^A^+yU + a). 
The assertion b) follows from a simple algebraic calculation. □ 



16 



REMARK 9 Since the Hilbert space is A^-dimensional with < 

oo, the resolvent formula ()28|) implies that the difference — G^^q, is a finite rank 
operator with rank less than or equal to A^. Thus the number, counting multiplicity, 
of negative eigenvalues of —A + e^A^ + is less than or equal to A^. 

Let us illustrate the approximation by point measures on a simple example in 
dimension two. Suppose that measure m is minus length measure supported by a 
circle of radius R, i.e. m is constant and negative measure. This makes the choice of 
approximating point measures very simple: we spread equidistantly points along 
the circle and all the points have the same coupling constant c 

727ri? 



Due to the symmetry, the spectrum of —A + m for this specific measure is known; 
it consists of the essential spectrum [0, oo) and a finite number of negative eigenvalues, 
which are all except the lowest one twice degenerate, see j2]. To find the eigenvalues, 
one has to decompose L^(M^) into angular momentum subspaces and then to look for 
solutions of an implicit equation in each of the subspaces. Therefore we can compute 
and compare both exact and approximate eigenvalues. 
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(a) (b) 

Figure 1: The dependence of the approximate eigenvalues on the number of point 
potentials for circle with R = 10 and e = 0.1 (a), e = 0.01 (b). The dashed lines 
represent the exact eigenvalues of —A + m. 
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Figure 2: The dependence of the approximate eigenvalues on the number of point 
potentials for R = 10, using the standard two-dimensional point potentials. The 
dashed lines represent the exact eigenvalues of —A + m. 

Each approximation is characterized by a pair of numbers, e > and N & N. In 
numerical calculations we fix e and we let N grow. The results for one chosen radius 
and two different parameters e are depicted in figurelH cases (a) and (b) correspond to 
e = 0.1 and e = 0.01, respectively. We observe that below some threshold number of 
points, the approximate discrete spectrum has no resemblance to the exact spectrum. 
The approximate eigenvalues may be very large negative and their number may be 
much higher than the number of exact eigenvalues (in figure ^ we even have not 
plotted all the eigenvalues which exist only for small N.) 

It appears that for larger e, we get a fast convergence of eigenvalues, however, 
they are all shifted from the exact ones. The reason is that since we work with fixed 
e, the limit operator is in fact —A + e^A^ + m instead of —A + m. On the contrary, 
small e means that one need more points to obtain a qualitatively correct spectrum, 
but then for a large number of points one gets much closer to the exact spectrum. 

We can also compare this approximation to [8', where approximating operators 
were Laplacians with point potentials. Those point potentials are of course different, 
they are not defined via a quadratic form and cannot be understood as a special case 
e = of section ini instead boundary conditons on wavefunctions are used, see p. 
Figure El presents the eigenvalues of Laplacians perturbed by point potentials which 
converge to — A + m with the same measure m as above. We have already mentioned 
in the introduction that here, we obtain a stronger convergence result than the one 
in jH]. Moreover, comparing both figures ^ and we see that employing fourth-order 
differential operators in the approximation may improve significantly the spectral 
convergence. 
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Appendix 

In section In] we have employed Krein's formula Various forms of this formula 
can be found in the literature. Let us prove here the one we have used. 

Let / G L^{dx). Since £e and are associated to —A + e^A"^ and —A + e^A'^ + fi, 
respectively, it follows from Kato's representation theorem that 

£eAGe,af, h) = (/, h) = S^Ji'A + A' + /i + a)-'f, h), (34) 

for any h E H'^iW'-) and / G L^(R°'). Moreover we have 

SUG^^.^, h) = {G^^ij, (-A + e'A' + a)h) 
= If QeA^ - y)i^iy)f^idy){-A + e'^A^ + a)h{x) dx 
= J J 9eAx - y)i-^ + ^^^^ + <y)h{x) dxtij{y)fi{dy) 
= Jhtpi^idy), ^eH^R"^), h e D{-A + e^A^). (35) 

We could change the order of integration in the second step. In fact, as /i^ are finite 
Radon measures and g^^a is square integrable w.r.t. the Lebesgue measure dx, the 
mappings x J [(/^^^(x — ?/)|/i^((i?/), M.^ — >• M, are square integrable w.r.t. dx. Since 
ijj is bounded and (—A + e^A^ + a)h G L'^{dx) it follows that 

9s A^ - y)i>{y)\iJ,^{dy)\{-A + e^A^ + a)h{x)\dx < oo 

and, by Fubini's theorem, we could change the order of integration in the second 
step. In the last step we have used Employing Sobolev's inequality and the 
fact that D^—A + e^A^) is dense in {D{£^),£^Ai we can extend to all functions 
ip,he D{£e). 
Put 

Then G H'^(R'^) = D{S^) and (jSl and extended ^ yield 
£^A^, h) = SUGs,af, h) - SeAGUl + G^,A''G,,J, h) 
+ j [G,,J - G^A^ + G^A''Gs,J]hd^i 

= {f,h)- l[{I + G^A''Gs,af]hdfi 

+ J[{i + G^ii + g':a-'g,,j - g>:ai + G^,A''Gs,J]hd^^ 
= if,h), heH'iR'). 

Due to (Uni), S^^^ is a scalar product on D{S^A = H^(^'^)- Thus §^ and the 
calculation above imply that = (—A + e^A^ + n + a)^^f. 




19 



Acknowledgment 

This work is partially supported by the Marie Curie fellowship MEIF-CT-2004- 
009256. 

References 

[I] S. Albeverio, F. Gesztesy, R. H0egh-Krohn, H. Holden: Solvable models in quan- 
tum mechanics, second edition. AMS Chelsea Publ. 2005. 

[2] J.-P. Antoine, F. Gesztesy, J. Shabani: Exactly solvable models of sphere interac- 
tions in quantum mechanics J. Phys. A 20 (1987), 3687-3712. 

[3] J. F. Brasche: On the spectral properties of singular perturbed operators, pp. 65- 
72 in Z. Ma, M. Rockner, A. Yan (eds.): Stochastic Processes and Dirichlet forms, 
de Gruyter 1995. 

[4] J. F. Brasche: Upper bounds for Neumann - Schatten norms. Potential Analysis 
14 (2001), 175 - 205. 

[5] J. F. Brasche, P. Exner, Y. Kuperin, P. Seba: Schrodinger operators with singular 
interactions. Journ. Math. Anal. Appl. 183 (1994), 112-139. 

[6] J. F. Brasche, R. Figari, A. Tcta: Singular Schrodinger operators as limits of point 
interaction Hamiltonians. Potential Analysis 8, no. 2 (1998), 163 - 178. 

[7] H. L. Cycon, R. G. Froese, W. Kirsch, B. Simon: Schrodinger Operators. Springer, 
Berlin-Heidelberg-New York 1987. 

[8] P. Exner, K. Nemcova: Leaky quantum graphs: approximations by point inter- 
action Hamiltonians. J. Phys. A 36 (2003), 10173-10193. 

[9] K. S. Murigi: On Eigenvalues of Schrodinger Operators with 5-Potentials. Master's 
Thesis, Mathematics, Chalmers University of Technology, Goteborg 2004. 

[10] A. Posilicano: Convergence of distorted Brownian motions and singular Hamil- 
tonians. Potential Analysis 5 (1996), 241-271. 

[II] P. StoUmann, J. Voigt: Perturbation of Dirichlet forms by measures. Potential 
Analysis 5 (1996), 109-138. 



20 



